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I/"") . Abstract 

oo ■ 

. We use the available information on the ratio between the electric and magnetic proton form 

£ — ' factors coming from recently published space-like data and from the few available time-like data. 

(~^) We apply a dispersive procedure on these data to evaluate the behaviour of this ratio, as a complex 

■ function, for all values of q 2 . 

o 

p ■ 1 Introduction 

Qh. The electromagnetic form factors (ff) are essential pieces of our knowledge of the internal structure 
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of the nucleon, and this fact justifies the efforts devoted to their determination. They are complex 
functions of the squared momentum transfer in the photon-nucleon vertex, defined both for space-like 
(q 2 < 0) and time-like (q 2 > 0) momenta. The interest in the study of the nucleon ff's has been 
recently renewed [1,2] by the Jlab unexpected results on the ratio 



R{q) -^Gwr () 

where [i v is the proton magnetic moment and G E (q 2 ) and G p M (q 2 ) are the electric and magnetic 
proton ff's. A decrease of R(q 2 ) as space-like \q 2 \ increases [3] has been found, in contrast with the 
flat behaviour R(q 2 ) ~ 1 [4] obtained up to q 2 f=a —7 GeV 2 by using the traditional Rosembluth 
method [5]. It has to be mentioned that the IJL nucleon ff's model [6], somehow related to soliton 
models of the nucleon [7], predicted the decrease of R{q 2 ) more than thirty years ago. Recently this 
model [8] and other QCD based models [1,9,10] have been extrapolated to time-like positive q 2 . While 
they agree within the experimental errors in the space-like region, they disagree in some respect for 
time-like q 2 . In all these models the space- like q 2 variation is related to a cancellation between the 
Dirac, Ff(q 2 ), and the Pauli, F^q 2 ), ff's. Then, since we have: 

G p E (q 2 ) = Ff(q)+rF!(q 2 ) 

T= 4^' (2) 
Gl I (q 2 )^Ff(q)+F p (q 2 ) AM P 

such a cancellation should become an enhancement once q 2 (i.e. r) changes sign from negative to 
positive, time-like values. As a consequence the angular distribution of the outgoing nucleon in 
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e + e~ — ► NN processes should have a large sin 2 (0) term proportional to \R(q 2 )\ 2 . At the same time, 
a large transverse polarization of the outgoing nucleon [11] is predicted [1], but the sign and the 
g 2 -dependcncc strongly depends on the models. 

In principle, time-like ff's could be evaluated from the space-like ones by means of dispersion 
relations (DR) , if they are smooth enough and if space-like data were known with very high accuracy or 
if in the time-like region there are some data or suitable constraints [12-14]. In fact, analytic functions 
are supposed to describe at the same time space-like and time-like electric and magnetic ff's. These 
functions are defined in the whole q 2 complex plane; in the physical sheet they do not have isolated 
singularities, such as poles, but a cut on the real axis starting at sq = (2m 1T ) 2 . Cauchy theorem 
allows to relate the space-like real values of a ff to an integral, over the time-like cut, of its imaginary 
part, providing a DR among them. In the unphysical region, < q 2 < &M 2 , not directly accessible 
experimentally, each ff should have large bumps corresponding to the unphysical excitations of p, u>, 
p', <J and all vector mesons whose mass is lower than 2M p [15]. In terms of ff's, these resonances are 
described by poles lying in the unphysical sheet of the q 2 complex plane. In the ratio R(q 2 ) the effects 
of these poles should be somewhat smoothed out, being the same in both, G E and G V M . Hence R(q 2 ) 
smoothness should be a plausible ansatz. 
There are constraints to R{q 2 ) for time- like q 2 , namely: 

• a continuous transition is expected from space-like to time-like q 2 ; 

• according to the Phragmen-Lindcloff theorem the space-like and time-like values of a ff must be 
asymptotically equal in modulus [16-18]; 

• at q 2 = si = 4M 2 , the physical threshold of e + e~ — > pp, assuming Ff(q 2 ) and F^{q 2 ) analytic 
functions of q 2 , it is R{s\) = fj, p . Since in this q 2 region R(q 2 ) is a complex function, this is a 
constraint on both, real and imaginary part. 

All the published data in the time- like region assumed |i?(g 2 )| = \x v not only at threshold 
[19-25], but in the whole explored interval [26-28], essentially for lack of accurate data concerning 
the angular distributions and of any measurement of the outgoing nucleon polarization or any use of 
polarized beams. 

In this paper, the ratio R(q 2 ) has been obtained in a model independent way, assuming 
it has a smooth behaviour, taking into account the theoretical constraints and solving the DR by a 
minimization algorithm. The input of this algorithm is given by the JLab and MIT-Bates polarization 
data [3,29] in the space-like region and by the result of a reanalysis of FENICE [19,20], DM2 [21,22] 
and E835 [30] data in the time-like region. A preliminary version [31] of this paper has been previously 
presented. 

A very similar approach, even if with a different dispersive integral, has already been success- 
fully tested on the pion ff [13]. In short, it has been assumed the pion ff to be known in the space-like 
region and in the time-like region, above Si only. Hence the DR has been reversed evaluating the 
time-like pion ff below si by means of a standard regularization procedure used to solve first kind 
integral equations, that requires one free parameter only [32]. A very good agreement has been ob- 
tained in this way, varying this parameter within one order of magnitude. This procedure has also 
been applied quite successfully to get the magnetic nucleon ff's in the unphysical region [13]. 

2 Space-like and time-like experimental data 

The high intensity-high polarization electron beams available at JLab and MIT-Bates [3, 29] allowed 
the extraction of the ratio between electric and magnetic proton ff's measuring the electron-to-proton 
polarization transfer. These measurements showed an almost linear decrease of R(Q 2 ) from unity 
at low Q 2 = —q 2 up to ss 0.3 at the highest Q 2 , as shown in Fig^ in strong disagreement with 
previous Rosembluth measurements (see rcf. [4] and ref. [33] for a compilation of all the space-like 
data), that indicated approximate ff scaling, i.e. R(Q 2 ) ~ 1, though with large uncertainties in G E 
at the highest Q 2 values. In the Rosembluth measurements, the ff's arc basically extracted fitting 
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the linear e-dependence of the cross section at fixed Q 2 (where e is the virtual photon polarization), 
but the presence of the factor 1/r [sec cq.©] in front of G P E (Q 2 ) makes difficult its extraction when 
Q 2 becomes large. Then, possible experimental inconsistency could affect the Rosembluth result, 
however a recent rcanalysis [33] showed that the individual Rosembluth measurements are consistent 
to each other within small normalization uncertainty between the different experiments. In addition, 
the new Rosembluth measurements performed at JLab [34] confirmed the scaling behaviour of the old 
data for Q 2 from 2.6 to 4.1 GeV 2 , making it clear that the source of the discrepancy is not simply 
experimental. 
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Figure 1: The full circles represent the data from Jlab and MIT-Bates [3,29], while the two time- 
like intervals are from FENICE [19,20], DM2 [21,22] and E835 [30]. The stars are the theoretical 
constraints and the shadowed area is the unphysical region. The label R(q 2 ) of the ordinate axis is 
defined as: R(q 2 ) = R(q 2 ) for q 2 < s and R(q 2 ) = \R{q 2 )\ for q 2 > s . 

Recent theoretical works [35] have suggested that terms, related to two-photon exchange corrections 
to the lowest order QED diagram, may result in incorrect determination of the ff's from the measured 
cross section, while polarization measurements should be less sensitive to such corrections. As shown 
in ref. [36], the Rosembluth data with the inclusion of two-photon exchange contribution agree well 
with polarization data for Q 2 between 2 and 3 GeV 2 , while there is at least partial reconciliation for 
higher Q 2 . Once reliable calculations of these corrections will be available, then polarization and cross 
section data can be consistently combined to extract the ff's without ambiguity. 

In our following analysis, we will make use of the ratio R{Q 2 ) as obtained from polarization 
measurements, since, as discussed, these data are less sensitive to higher order corrections and to 
systematic uncertainties. 

The proton ff's in the time-like region (q 2 = s > 0) can be extracted from the cross-section of 
the process e + e~ — > pp or pp — > e + e~. The measurements available from both kinds of experiments 
are shown in Fig[21as a function of the center of mass energy s. Most of the data are concentrated in 
the low s region close to the proton-antiproton threshold [19-25], but a sizeable amount of data from 
proton-antiproton annihilation experiments is also available for center of mass energies s between 8 
and 14 GeV 2 [26-28]. The total cross-section at a given s is related to the ff's | (s) | and | (s) | 
through the relation: 
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where the factor p(s) is equal to C/3 = Cv/c (v is the proton velocity in the center of mass system) 
in the case of e + e~ — > pp and to f3/C in the case of pp — > e + e~ and C is a Coulomb correction factor 
to take into account QED bound states, relevant only very near threshold. In order to extract the 
ff's from the measured cross-sections, each experiment has to make an hypothesis on the modulus of 
the ratio R(s). All the results published are obtained in the hypothesis that this ratio, in modulus, is 
equal to Mp (\G p E (s)\ = \G p M (s)\). 
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Figure 2: Proton magnetic ff in the time-like region extracted from the e + e~ — ► pp and pp — > e + e~ 
cross-sections assuming |G^(s)| = |G^ f (s)|. a) All data including the largest s ones in logarithmic 
scale, b) Data close to si = 4M 2 in linear scale. 



The possibility to disentangle between |G^(s)| and |G^(s)| relies on the measurement of the angular 
distributions. Calling 9 the polar angle of the emerging proton (or antiproton) in the center of mass 
system in e + e~ — * pp experiments, or the polar angle of the electron (or positron) in pp — > e + e~ 
experiments, the differential cross-section is 

% = ^ (Vm(*)| 2 (1 +cos 2 (fl)) + V £ ( S )| 2 sin 2 (0)) , (4) 

where r = s/4M 2 . This formula shows that the two ff's give rise to two terms: one, related to the 
magnetic ff, has a [1 + cos 2 (6>)] dependence, the other one, related to the electric ff has a sin 2 (6>) 
dependence. In principle a measurement of the angular distribution is able to give the relative weights 
of the two terms and hence the ratio |-R(s)|. A fit of the cos(0) distributions has been applied to two 
sets of data corresponding to two different values of s. 

1. The data from the FENICE [20] and DM2 [22] experiments both detecting the process e + e~— >pp 
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at an average center of mass energy s = 4.51 GeV 2 have been simultaneously fitted. The two 
measured differential cross-sections are in good agreement as shown in Fig|3Ji. 

2. Two sets of data from the E835 experiment [27, 28] detecting the process pp — > e + e~ both at 
s ~ 12 GeV 2 have also been fitted simultaneously. The angular distributions of the two data 
sets are shown in FigEt normalized to the average number of events per bin. 
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Figure 3: a) Differential cross-section as a function of |cos(0)| from FENICE [19,20] (triangles) and 
DM2 [21,22] (full squares) data at s^ xp = 4.5 GeV 2 , with fit superimposed, b) Angular distributions for 
E835 1996 data (full squares) and 2000 data [30] (triangles) at s^ p = 12.0 GeV 2 , with fit superimposed: 
the two distributions are normalized to the same number of events per unit of cos(0). 



All the considered experiments are characterized by limited statistics due to the small value of the 
cross-section (from ~ 1 nb at s^ xp = 4.5 GeV 2 down to ~ 1 pb at s^* p = 12 GeV 2 ). In both cases the 
parameter |i?|^ xp (k = 1, 2) is extracted from a fit to the cm. angular distribution with: 

A(cos0) = rn 2 p A 2 k (l + cos 2 0) + B 2 k {l - cos 2 9), (5) 

the ratio between the two coefficients B k and A k is just |i?|^ xp - In order to evaluate confidence 
intervals for |-R|)* P properly taking into account the unavoidable correlations, the same fits have been 
applied on samples of distributions randomly extracted from the experimental ones. This procedure 
allows to get directly the distributions of |-R|^, xp - These are not gaussian, being in both cases strongly 
asymmetric. A 68% confidence intervals (-Rl < \R\ e k P < R 3 k up ) are built in such a way that the 
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probability P(\R\ > R s k up ) = P(\R\ < R^ } ) = 16%. The resulting intervals are: 

3.46 < \R\? P < 26.5 (sf p = 4.5 GeV 2 ) 

0.34 < |i?|2 xp < 3.63 (s^ = 12.0 GeV 2 ). W 

These numbers could be modified by taking into account the contribution of two-photon exchange 
diagrams, that should lead to an asymmetry in the angular distribution of the proton with respect to 
the antiproton. An effect of the order of few percent could be estimated, also taking into account that 
the 77 — ► pp cross section is of the same order of magnitude as for e + e~ — > pp [37]. However, given 
the experimental uncertainties in the determination of the two intervals of eq.©, we do not expect 
that our final result will be affected by these higher order corrections. 



3 The dispersive approach 

In order to connect the data on |i?(s)|, in the time-like region to those on R(t), in the space-like region 
(t = q 2 < 0), the DR's for the imaginary part is a very powerful mathematical tool that we can use. 
The DR for the imaginary part has the form [16]: 

It establishes an integral relation between the values of the ff in the space-like region, where this 
function is real since it describes the scattering process, and its imaginary part in the time- like region, 
where the process described is the annihilation and therefore the ff is complex. If we assume that 
the magnetic proton ff G p M (q 2 ) has no zeros, as it is demonstrated in ref. [13], it follows that, unless 
the asymptotic behaviour, the ratio R{q 2 ) has the same analytic properties of each ff. This means we 
may apply the DR of eq. (UJ) directly to R{q 2 )- But, since pQCD [17] constrains the ff 's G p E (q 2 ) and 
G p [I (q 2 ) to be asymptotically vanishing with the same power law (1/q 2 ) 2 , as q 2 diverges, the ratio, 
unless logarithmic corrections, should have a constant asymptotic time-like limit and therefore the 
integral of eq.Q), with R(s) instead of G(s), could be divergent. In order to account for this possibility 
we perform the analytic continuation of R(q 2 ) by means of the DR for the imaginary part, subtracted 
at t = [16]: 

Kj So S{s-t) 

which relates the space-like value of R(t) to its time-like imaginary part and: 

r , si , , s f°° lm\R(s')} , 
Rc[R( S )} = R(0) + -Vj^ p^MJd*', (9) 

which, instead, connects the real and the imaginary parts of R(s) over the cut, i.e., for s > sq (V 
denotes the principal value). The price of the subtraction at t = 0, is the knowledge of the value of 
R{t) at the same point, but, thanks to the normalization [eq.JTJ], in the expressions of the eas. l|8l9|l 
we can put i?(0) = 1. 

Contrary to the existing models [1,6,8-10], which are constructed starting from the space- like data 
and only subsequently extended to other energies, we start from the imaginary part of the ratio R(s), 
which is defined only in the portion of the time-like region over the cut, and, by means of a rigorous 
analytic continuation procedure, we reconstruct the function R(q 2 ) in the whole q 2 complex plane. 
To parametrize the imaginary part we use a quite general and model independent form, that is, two 
series of orthogonal Chcbyshcv polynomials [16]: 

ZfC&ix) x= 2 -^^ so<s<si 
lm[R(s)] = I(s) = { (10) 
EfDjT^y) y=^-l a > a u 
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with s > So, the two vectors C= (Ci, C2, . . . , Cm) and D= (D\,D2, ■ ■ ■ , Dn) represent the coefficients 
and Tj(x) is the j-th Chebyshev polynomial. The two series in ea. ljTU)) cover two naturally separated 
intervals: 

• the unphysical region [sq, s\] with the vector meson resonances, where the ff's can not be mea- 
sured; 

• the experimentally accessible region (si,oo), where the asymptotic regime is attained. 




q 2 (GeV 2 ) 

Figure 4: a) Result of our dispersive technique for the ratio R(q 2 ) in the full data region. The solid 
band represents the error, b) Ratio R{q 2 ) in the time-like region with q 2 > si and comparison among 
various models (Q 2 = —q 2 )- As in Fig^thc full circles arc the data from Jlab and MIT-Bates, while 
the two time-like intervals are from FENICE, DM2 and E835. The stars represent the theoretical 
constraints. The label R(q 2 ) of the ordinate axis is defined as: R(q 2 ) = R(q 2 ) for q 2 < so and 
R{q 2 ) = \R(q 2 )\ for q 2 > s . 
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Figure 5: a) Ratio G p E (q 2 ) / G p M {q 2 ) over a wide range of q 2 . The solid band represents the error, 
b) Our result compared with the other models (same references of Fig^J) in logarithmic scale for q 2 . 
The label G p E (q 2 )/Glj(q 2 ) of the ordinate axis is defined as: G p E {q 2 )/G p M (q 2 ) = G p E (q 2 ) / G p M {q 2 ) for 
q 2 < so and G p E (q 2 )/G p M (q 2 ) = \G p E (q 2 ) /G p M (q 2 )\ for q 2 > s . 

Once Im[i?(s)] is defined, R(q 2 ) can be evaluated in the full q 2 range via the eas. (|8l9|) . In order to find 
the C and D vectors we minimize a \ 2 function that includes space-like and time-like data together 
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with the constraints discussed above. The resulting function has the correct analytic structure. More 
details on the minimization procedure are reported in rcf. [38] . 

Concerning the dispersive procedure, it is interesting to mention that, in addition to the one 
adopted here, there are other forms of DR's which connect space-like and time-like data. For instance, 
the DR for the logarithm may relate directly data on the modulus of a ff in the time-like region to 
those for its real value in the space-like region [12,13,16]. This kind of logarithmic dispersive approach 
would not be suitable in this case, not only because the time-like modulus has not so many constraints 
as the imaginary part does, but also because it is not able to forecast, in a model independent way, 
the presence of a zero, that, instead, should be one of the aims of this analysis. In the pictures of 
Figgis shown our result for the ratio R{q 2 ). We found, in a model independent way, a space-like zero 
at q 2 = (—11 ± 2) GeV 2 (FigQJi). The comparison between some existing models [1,6,8-10] and our 
result, reported in FigQJ), shows that, in spite of the agreement in the space- like region, where all these 
models describe the polarization data, their continuations in the time-like region are far each other. 
Another interesting outcome of our computation is the space-like and time-like asymptotic behaviour 
of G E (q 2 )/G p M (q 2 ). According to the Phragmcn-Lindeloff theorem, that we have implemented by 
imposing a vanishing asymptotic value for the imaginary part of R(q 2 ), we achieved a real time-like 
limit for the ratio as q 2 — > +oo. As shown in Fig[S] (without and with comparison with the other 
models) the space-like and time-like asymptotic limits have, in modulus the same value, but they have 
opposite sign, i.e.: 



lim 



±1. 



(11) 



The scaling law, providing G P E ~ G p Ml that, in light of the new polarization data, is no more valid at 
low |g 2 |, has been in some way restored, even if in modulus and in a "double" asymptotic regime. In 
addition, the Phragmcn-Lindeloff theorem predicts for the phase 

$(oo) = [Z - P)tt, (12) 

where Z ans P are the number of zeros and poles of R(q 2 ). Since we have: R(—oo) = |i?(oo)|e** ( ' 00 ) = 
— |i?(oo)|, and being P = 0, we obtain a further confirmation of the presence of a zero of R(q 2 ) (or at 
least an odd number of zeros). 
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Figure 6: Ratio \F p (q 2 ) / Ff (l 2 )\, our prediction (solid band) compared with the other models (same 
references of FigQJ. 



As a direct consequence of the limit of ea. l|ll|l . we give also a prediction (FigjSJ) for the ratio between 
Pauli and Dirac ff 's as: 



T F jS q l = -0.2 ±0.3, 



— >oo 



Ff(q 2 



(13) 
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where, as usual, r = jfyy, therefore Ff \q 2 ) / Ff (q 2 ) scales at least like (1/g 2 ) as q 2 diverges. 

p 

The complete knowledge of the function R(q 2 ) allows to give predictions concerning also the quantities 
depending explicitly on real and imaginary part (or on modulus and phase) as polarization observ- 
ables. In the space-like region, the scattering of polarized leptons on a proton target gives non trivial 
polarization effects even if the target is unpolarizcd. The polarization of the outgoing proton, in this 
case, depends on the product of electric and magnetic ff 's, which in this region are real. 
In the time-like region, i.e. when we consider the annihilation e + e~ — > pp, the complex structure of 
the ff's give rise to special polarization effects: the outgoing proton may experience a polarization 
even if there are no polarized leptons in the initial state. This polarization will be determined by the 
relative phase of G p E {q 2 ) and G p M (q 2 ), that, in our case, is just the phase <E>(g 2 ) of R{q 2 )- The time-like 
polarization vector V has the components [11]: 



sin(260|i?(< ? 2 )|sin($(<7 2 )) 



Vy{q 2 ) 

2sin(6»)|i?(q 2 )|cos($(g 2 )) 



V x (q 2 ) = -p e ^^l 



2cos(6») 



D 1 
with: 

l + cos 2 (60 + ±|T| 2 sin 2 (<9) , N 

D= y — \1 } (15) 

where z is the direction of the outgoing proton in the center of mass system and y is orthogonal to the 
scattering plane, P e is the longitudinal polarization of the initial lepton and 9 is the scattering angle. 
As already said, the y-polarization V y does not depend on P e , while the longitudinal polarization 
does not depend on the pahse <&. In the pictures of Fig[7|are shown our predictions and those of the 
other considered models [1,6,8-10], for the components of the polarization vector. Once again models, 
which agree in the space-like region, give different predictions in the time-like region. 




q 2 (GeV 2 ) q 2 (GeV 2 ) q 2 (GeV 2 ) 

Figure 7: Our predictions (solid band) for polarizations V x , V v and V z in the time- like region compared 
with those of the other models (same references of FigQJ . The plots are for the scattering angle 8 = 45° 
and electron polarization P e = 1. 



4 Conclusions and perspectives 

We have used a dispersive approach to construct an expression for the ratio R(q 2 ), defined in the whole 
q 2 complex plane, which verifies all the constraints imposed by the theory (analyticity, asymptotic 
behaviour, etc.) and by the data available at this moment. The full knowledge of the function R(q 2 ) 
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allows to formulate a wide range of predictions. In particular, in the time like region, dominance of 
the electric form factor near threshold, fading away soon, as well as an oscillating pattern is predicted 
(see Fig0J, that could be interpreted as a resonances residual, survived to the smoothing expected 
in the ratio. Also the ratio between Dirac and Pauli ff's is predicted as well as a definite nucleon 
polarization. The other prediction that should be confirmed or refuted soon, is the presence of the 
space-like zero, that we estimate at q 2 = ( — 1 1 ± 2) GeV 2 , in agreement with ref. [39]. New polarization 
measurements are scheduled at Jlab to push down the space- like limit at q 2 ~ — 10 GeV 2 [40]. 
Anyway, the measurement of |i?(q 2 )| and of the polarizations [eq. H14|l ] in the time-like region would 
have a crucial significance not only to disentangle among the models, but also to gain a rather complete 
experimental knowledge of the proton ff's. 
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